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Abstract. Modulo conventional scale factors, the Simon and Simon-Mars 
tensors are defined for stationary vacuum spacetimes so that their equality follows 
from the Bianchi identities of the second kind. In the nonvacuum case one 
can absorb additional source terms into a redefinition of the Simon tensor so 
that this equality is maintained. Among the electrovacuum class of solutions 
of the Einstein-Maxwell equations, the expression for the Simon tensor in the 
Kerr-Newman-Taub-NUT spacetime in terms of the Ernst potential is formally 
the same as in the vacuum case (modulo a scale factor), and its vanishing 
guarantees the simultaneous alignment of the principal null directions of the Weyl 
tensor, the Papapetrou field associated with the timelikc Killing vector field, the 
electromagnetic field of the spacetime and even the Killing-Yano tensor. 



PACS number: 04.20.Cv 



1. Introduction 



In the 1 + 3 gravitoelectromagnetic approach to general relativity associated with 
splitting spacetime quantities using a unit timelike vector field u, the 4-velocity of an 
observer congruence, the Bianchi identities for the curvature tensor take a Maxwell- 
like form for the electric and magnetic parts of the Weyl tensor with sources arising 
from the Ricci tensor terms, which in turn can be expressed in terms of the energy- 
momentum tensor of the spacetime through the Einstein field equations [1-3]. For 
the stationary case with u directed along the timelike Killing vector field, these 
equations are greatly simplified and suggested the introduction of the complex Simon 
tensor [4], which in the vacuum case equals the simpler Simon-Mars tensor [5-10] 
(modulo conventional scale factors). The Simon-Mars tensor is a natural algebraic 
combination of the self-dual fields associated with the Papapetrou field of the timelike 
Killing vector field and the Weyl tensor. 
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Simon introduced his Simon tensor for vacuum stationary spacetimes as a complex 
generalization of the Cotton- York tensor, to which it reduces in the static case, where 
its vanishing guarantees that the space sections are conformally flat, as occurs in the 
Schwarzschild spacetime. Its vanishing also characterizes the Kerr spacetime within 
certain classes of solutions of the stationary vacuum Einstein equations. Mars has 
given an equivalent form of the Simon tensor (the Simon-Mars tensor) that has been 
very useful in understanding the meaning of its vanishing in the Kerr spacetime, 
where it leads to the alignment of the principal null directions of the Weyl tensor of 
the spacetime with those of the Papapetrou field associated with the timelike Killing 
congruence. 

Here we extend the gravitoelectromagnetic analysis of the vacuum case [10] to 
electrovac stationary spacetimes which contain a source-free electromagnetic field, 
where the Simon tensor and the Simon-Mars tensor defined by the same formulas 
as in vacuum differ by additional source terms. However, these source terms can be 
absorbed into a redefinition of the Simon tensor so that it continues to be equivalent 
to the Simon-Mars tensor defined exactly as in the vacuum case, and of course having 
all the same tensor symmetries. Remarkably for the Kerr-Newman spacetime and 
even the more general family of Kerr-Newman- Taub-NUT spacetimes, the situation 
is very similar to the Kerr spacetime. First the vanishing of the Simon-Mars tensor 
(and hence of the newly defined Simon tensor) can still be explained in terms of 
the simultaneous alignment of the principal null directions of the Weyl tensor, the 
Papapetrou field associated with the timelike Killing vector field, the electromagnetic 
field of the spacetime and even the Killing- Yano tensor. Secondly when written in 
terms of the Ernst potential, the newly defined Simon tensor including the source 
terms formally coincides with the vacuum expression (modulo a scale factor). 

Our extension of the Simon tensor to stationary electrovac spacetimes takes a 
different approach from the work of Marklund and Perjes on stationary perfect fluid 
spacetimes [11]. They are also able to incorporate additional source terms into the 
definition of the nonvacuum generalization of the Simon tensor in such a way that 
the symmetry properties of the tensor are respected. However, they focus on the 
spatial Bianchi identities (for the conformally rescaled spatial metric) rather than the 
spacetime Bianchi identities. 

Finally it is worth noting that a prime motivation of our discussion, begun 
in [10], is to make more accessible the geometry surrounding the Simon tensor, 
whose existing literature assumes detailed knowledge of stationary exact solution 
techniques that outsiders simply do not possess. Nonexperts cannot even recognize 
the equivalence of different formulas for the Simon tensor stated by various experts 
in these manipulations. Appendix A reconciles these definitions, while the main body 
recasts the geometry in the natural complex version of the Maxwell-like equations for 
the Weyl tensor pioneered by Ellis [1-3]. 

2. The 1 + 3 Maxwell-like equations for the Weyl tensor with sources 

The Bianchi identities of the second kind equate the divergence of the Weyl curvature 
tensor C a /3 7 s to half the Cotton tensor [12], which can be thought of as a current in 
analogy with electromagnetism [13] (where F a/3 = A-kJ") 

C aS fj-,-8 = -V[ (3 (i? Q 7 ] - -i?(5 Q 7 ]) = J a (il , (2.1) 
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following the conventions of Misner, Thorne and Wheeler [14], including the signature 
— h++. The Einstein equations (in Ricci form) 

R a p = <T a - \T\5 a p ) (2.2) 

can then be used to replace the Ricci tensor and scalar curvature terms in this current 
by the energy-momentum tensor. Splitting these equations with respect to a generic 
timelike congruence u leads to their 1 + 3 Maxwell-like form [3] given explicitly in our 
notation in [10]. 

The Weyl tensor splits into two symmetric tracefree spatial fields, its electric and 
magnetic parts respectively 

E a = C a ^s^u s , H a = -*C a 70S uiu 5 = ^ V a 7 5 C">sp P u p . (2.3) 

If a is the acceleration vector and u> the vorticity vector of u, introducing the complex 
spatial fields from the gravitoelectromagnetic connection and Weyl curvature fields 

z = -a-iuj = g-iH/2 , Z = E — iH , (2.4) 

where g and H are the gravitoelectric and gravitomagnetic vector fields, leads to a more 
efficient representation of those equations. These complex fields are associated with 
the 0(3, C) representation of the Lorentz group exploited in the 'complex vector / self- 
dual formalism' [15-17] in turn closely related to the Newman-Penrose formalism. Here 
we specialize to the case of a stationary electrovacuum spacetime and align u = M^ 1 ^ 
with the associated Killing vector £, so that the acceleration can be expressed as 
a = VlnM, where M = | -£ Q £"| 1/2 . 

Splitting the Weyl divergence equations and expressing them in terms of an 
adapted frame ey = u,e a , a = 1,2,3, where the spatial frame e a is orthogonal to 
u, one finds in an index-free notation the identities [10] 

M 3 div[M- 3 Z] = 3z-Z + p^ , M- 1 Scurl[MZ] = z x Z + i J (G ^ , (2.5) 

where the complex current fields are defined by 

P ( a G) = J T aT + *J* T aT , = -[J (a b)T + iJ* (ab)T ] . (2.6) 

The div and Scurl spatial operators on symmetric spatial 2-tensors S 

[div(u)S} a = \I{u)pS af3 , [Scurl(ii)S*] Q/3 = n{u)~> s{a \7 {u)^ s (2.7) 

are defined in terms of the spatial covariant derivative V(w), the spatial volume 3- 
form 77(u) Q /3 7 = u^nsap-y (associated with the spatial duality operation *<">, used for 
example to define the vorticity vector u)(u) a = i?7(u) a/37 V / 3U 7 from the corresponding 
2-form, see appendix A of [10] for sign conventions) and the spatial cross and dot 
products associated with u, all given explicitly in [10] by spatially projecting spacetime 
derivatives and quantities orthogonally to u. For example, the cross-product of a 
spatial vector and spatial symmetric tensor is 

[X x u A] afj = ri(uy s ( a X^A® s . (2.8) 

For simplicity we drop the reference to u on these spatially projected operators, 
resorting to the abbreviation V = V(u) to distinguish the spatial covariant derivative 
from the spacetime covariant derivative V, an unnecessary distinction when acting 
on stationary scalars. In an adapted frame spatial quantities may be expressed using 
only Latin indexed components. 
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The spatial cross product of z and Z which appears in the Scurl equation in 
the projected Bianchi identity l|2.5[l defines the (complex, symmetric and tracefree) 
Simon-Mars 2-tensor 

SimonMars = z x Z . (2.9) 

By using that Scurl equation to define the nonvacuum Simon 2-tensor for any 
stationary spacetime (with any energy-momentum source) as 

Simon = M' 1 Scurl[AfZ] - i J (G) , (2.10) 

this identity then guarantees that the Simon and Simon-Mars tensors also remain 
equal in the nonvacuum case. Of course once the Ricci tensor in the source term j( G ' 
is replaced using the Einstein equations (|2.2|l , this equality between Simon and Simon- 
Mars is only true if they are each evaluated on a solution of those equations. Apart 
from a scale factor, our Simon tensor is the spatial dual of the three-index spacetime 
tensor introduced by Simon in the vacuum case [4], as shown in Appendix A. 

In the special case in which the source term has the form j( G ) = — iVlncr x Z, 
then the Scurl identity for any scalar a and any symmetric spatial 2-tensor A 

crScurl[<7 _1 vl] = -Vino- x A + Scurl A (2.11) 

enables it to be incorporated into the existing Scurl term 

Simon = (a/M) Scurl [cr^A/Z] . (2.12) 

This again would re-establish the vacuum property that the Simon tensor is 
proportional to a Scurl, but with a generalized Simon tensor incorporating the source 
terms so that it still equals the Simon-Mars tensor. This is the case with the family 
of Kerr-Newman-Taub-NUT electrovac spacetimes (including the more physically 
interesting Kerr-Newman spacetime), where the extra factor of a converts this to 
the same Scurl expression when written in terms of the Ernst potential as in the 
vacuum formula (modulo a scale factor). 

The vanishing of the Simon-Mars tensor (z x Z = 0) is equivalent to the alignment 
of z and Z in the sense that Z must be proportional to the tracefree tensor product 
of z with itself [10] 

Z = k[z® z]( TF ) , (2.13) 

which is an algebraic consequence of the definition of the cross product for symmetric 
spatial 2-tensors [10]. Thus the vanishing of the generalized Simon tensor corresponds 
to the alignment of z and Z as in the vacuum case. 

3. Maxwell equations and the Papapetrou field 

The stationary Maxwell-like equations i|2.5|) for the Weyl tensor expressed in complex 
self-dual form are very similar to the stationary Maxwell equations expressed in 
complex form [10] 

M~ 2 &\v{M 2 x) = -2z ■ x + 4tt P , M^ 1 cml(Mx) = iniJ (3.1) 

satisfied by the spacetime's electromagnetic field (using an arrow to distinguish the 
electric field from the electric part of the Weyl tensor in index-free notation) 

F = uAE + *^B, x = E-iB. (3.2) 

They only differ in form by the missing zxx term in the curl equation here, where *WB 
is the spatial dual of the magnetic 1-form field B. However, when z and Z are aligned 
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(i.e., z x Z = 0) as in the Kerr-Newman-Taub-NUT spacetimes, the correspondence 
is complete. The sources p and J for the spacetime electromagnetic field are zero for 
the electrovac spacetimes considered here. The curl equation for the electromagnetic 
field implies that Mx is the gradient of a potential 

Mx = V$ , or x = M -1 V$ . (3.3) 
Papapetrou [18] observed that a Killing vector field £ can be interpreted as a 
Lorentz gauge vector potential for an electromagnetic field later called the Papapetrou 
field [8,9] 

F a /3 = [dftap = V a £p - = 2V Q £/5 (3.4) 

satisfying Maxwell's equations with an effective current source proportional to the 
Ricci tensor 

= - = -> ^ = ~^BP a -» £,/3 i/? = 2&i^ Q , (3.5) 

where this divergence condition is merely the trace of the Ricci identity for a 
Killing vector field. Rescaling the Papapetrou field F by the factor (2M) _1 [10] 
and introducing its complex self-dual representation yields the self-dual 2-form T 
associated with z 

T = (2M)- 1 (F + i *F) = «Az + i*^z , (3.6) 

which when multiplied by M satisfies the same Maxwell equations but with nonzero 
sources generated by the electromagnetic field. These equations 

Af~ 3 div(M 3 z) = -2z ■ z - x ■ x , M" 2 curl(M 2 z) = x x x, . (3.7) 
are obtained by splitting the vector Eqs. (|3.5() and using the Einstein equations l|2.2|l 
to replace the Ricci tensor appearing there in terms of the appropriate projections of 
the electromagnetic energy-momentum tensor of the spacetime 

kT t j- = —x ■ x , inT T a = [x X x] a , (3-8) 

which then leads to the source terms appearing in these Maxwell equations. 

Inserting the potential representation (|3.3(l of Mx into the rescaled Papapetrou 
field curl equation of Eqs. Ij3.7|l then yields 

M~ 2 curl(M 2 z) = M" 2 V$ x V# = —M~ 2 curl($V$) , (3.9) 

namely 

curl(Af 2 z + $V$) = . (3.10) 

This latter quantity in parentheses must therefore be a gradient, defining the Ernst 
potential £ [15] (see their Eq. (18.34) with F = —M 2 , reducing to the integrated 
Eq. (18.29) in the vacuum case $ = 0) modulo a numerical factor 

V£ = -2(M 2 z + $V$) . (3.11) 
The remaining divergence equation of Eqs. I|3.7|l is the Ernst equation, which can 
be expressed in the form (equivalent to Eq. (18.38) of [15], see Appendix A) 

M div(MW) = -2M 2 z • V£ . (3.12) 

This result follows from using the definition of £ to re-express x in terms of the gradient 
of <!> in the two divergence equations (|3.1|) and (|3.7|) 

- div V£ + i>divV$ = M 2 (a + 2z) ■ z, div V$ = -{a + 2z) ■ V$ (3.13) 

and using the second of these equations to eliminate div in the first, yielding 

divV£ + a- VS = -2z- V£ , (3.14) 
which can then be rewritten as Eq. (|3.12() . 
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4. Simon and the Ernst potential 

The Simon-Mars tensor can be evaluated in a form which shows its relation to the Ernst 
potential expression for the Simon tensor, assuming that the Einstein field equations 
are satisfied so that the two tensors are in fact equal. Following the notation of [10], 
the purely spatial components and the mixed time-space components of the Ricci 
tensor field equations (see [15], Eqs. 18.15 and 18.16) are respectively 

[WiUF F) = -2M- 2 [V (a $V fc) $]( TF ) , 

^R Tb =iM- 2 [curl(M 2 z)] b = -iM- 2 [curl($V$)] fc , (4.1) 

where z and <& are related to the Ernst potential by (|3.11|) . A straightforward 
computation gives the electric and magnetic part of the Weyl tensor in terms of 
the gravitoelectromagnetic fields a and u> and the Ricci tensor of the quotient 3- 
manifold [10] 

E = - i[4w®w-Va-a®a- (3) Ricci] (TF) , 

H = - SYM[Vw + 2a ® w] (TF) . (4.2) 
Then using the formula 

[( 3 > Ricci] < TF ) = [^Ricci + Va + a®a + 2LU® w] (TF ) , (4.3) 

one finds 

Z = E -iH =[- SYM(Vz) + z <g> z + i (4) Ricci] (TF) . (4.4) 

The Simon-Mars tensor is then 

SimonMars = z x (- SYM Vz + i (4) Ricci) (4.5) 

since z x (z ® z) is identically zero. By next replacing the first term here using the 
Scurl identity valid for an arbitrary spatial vector field X 

Scur\(X ®I)=-Ix SYM(VX) + ^[SYM(X ® curl Jf )] (TF) , (4.6) 

and using the Einstein Eqs. I|4.1|) to replace the spatial components of the spacctimc 
Ricci tensor, the Simon-Mars tensor becomes 

SimonMars = Scurl(z ® z) - -[SYM(z ® curl z)] (TF) 

- M~ 2 z x SYM(V$ ® V<J>) . (4.7) 

Finally by substituting z — M~ 2 (M 2 z) into the curl term here expanding the 
result with the appropriate product rule and incorporating the resulting additional 
acceleration term into the Scurl term, and then using Eq. H3.11|l 

curl(Af 2 z) = - curl($V$) (4.8) 

to substitute the other term in this expansion, one finds 

SimonMars = Af~ 3 Scurl(M 3 z®z)- ^[SYM(z <g> (x x s)] (TF) 

- z x SYM (a; ® x) . (4.9) 
In the vacuum case where x — 0, this reduces to 

SimonMars = Af ~ 3 Scurl(M 3 z ® z) , (4.10) 
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which can be rewritten in terms of the Ernst potential by making the substitution 
z = -\M- 2 V£ 

SimonMars = Jm~ 3 Scur^M - ^ <8> V£) ■ (4.11) 

As shown in appendix A, this is proportional to the form of the Simon tensor given 
in [15] and to the original form given by Simon. 

5. The Kerr-Newman-Taub-NUT spacetime 

In Boyer-Lindquist coordinates (x a ) — (t,r,9,(p), with the abbreviations (c, s) = 
(cos 9, sin 9), the (exterior) Kerr-Newman-Taub-NUT (KNTN) spacetime the line 
element is [19-23] 

d.s 2 = - i(A - aV)di 2 + | [A X - a(S + a X )s 2 }dtd<j) 

+ i[(S + a X fs 2 - X 2 A]d0 2 + |dr 2 + Zd9 2 , (5.1) 

and the corresponding electromagnetic Faraday tensor can be expressed in terms of 
the 2-form 

F = ^{[r 2 -{£ + ac) 2 }dr A (dt - X d0) 

+ 2rs(£ + ac)d9 A [{r 2 + a 2 + £ 2 )d<j) - adt]} , (5.2) 
which is the exterior derivative of the vector potential 

A = -^-(dt- x d<j>) . (5.3) 

Here S, A, and \ are defined by 

S = r 2 + (£ + ac) 2 , A = r 2 - 2Mr - i 2 + a 2 + Q 2 , X = as 2 - 2£c . (5.4) 

Units are chosen so that G = c = 1, so the parameters (M,Q,a,£) all have the 
dimension of length. The source has mass M, electric charge Q, angular momentum 
J = Ma (i.e., gravitomagnetic dipolc moment) along the 9 = direction, and 
gravitomagnetic monopole moment fj, = —£, where £ is the NUT parameter. [24-27]. 

The spatial 3-metric h a p on the quotient space of the timelikc Killing congruence 
associated with the Killing vector field £ = d t is 

d^ 2 = h ap dx a dx fi = ^dr 2 + £d# 2 + A As2 ^ 9 d0 2 , (5.5) 
A A — a z s z 

while its convenient conformal rescaling by the square of the lapse function M 2 = 
= (A - a 2 s 2 6»)/E is 

df 2 = ~h a[j dx a dxP = M 2 h a[i dx a dx' 3 

= A -^ Sin2g dr 2 + (A _ a 2 s 2 )d() 2 + As 2 d(f> 2 (5 g) 

Using the latter metric and its covariant derivative V, one easily shows by direct 
evaluation that both the Ernst potential and $ 

g = l 2 M+3 * = g (57) 
r + i(£ + ac) ' r + i{£ + ac) 
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satisfy the same Scurlfree property for their gradient 

Scurl (V£ <E> V£ ) = (5.8) 

that occurs for the Ernst potential in the vacuum case within this family. 
This is explained by the special form of the Weyl current 

J(G) = [r - i(l + ac)] [r + i{£ + ac)] 2 W 2 (A - a 2 s 2 fl 2 ^ ~ A ^ 

= -iVloger x Z , (5.9) 

where 

3 {[r - i(£ + ac)](M + it) - Q 2 } [r + i(l + ac)} 
°~ 4 {M+iif [r-i(e + ac)} ( ' 

a factor also which appears in 

Z = M- 2 a[V£ <g> V£ ] (TF) , (5.11) 

and so cancels the factor inside the Scurl in the Simon tensor expression H2.12|l when 
expressed in terms of £ 

Simon = crM" 1 Scurl(cr _x MZ) = aM' 1 Scurl(M _1 VS ® V£) . (5.12) 

This latter Scurl expression is the same one originally used to define the Simon tensor 
in the vacuum case (where Vcr x Z = and so disappears from the formula upon 
expansion of the first Scurl expression), and as shown in Appendix A is proportional 
to 5curl(V£ (?) V£). This is forced to be zero when z and Z are aligned as they are 
here 

{[r - i{t + ac)]{M + it) - Q 2 } 



[r + i(l + ac)] 2 (A - a 2 s 2 )[r - i(l + ac)] 



{-Ad r + iasd e ) , (5.13) 



Z - 3(A-aV) (TF) 

In fact all the complex vectors z, V£, V$, and the one associated with the Killing- Yano 
tensor are all aligned with (i.e., proportional to) the gradient vector field W[r+i(£+ac)]. 
Although we have downplayed the notational distinction between contravariant and 
covariant tensors, by writing explicitly the index-lowered 1-form 

Z (A- a 2, s 2) [r + i(l + ac)] 2 [r _ i{l + ac)] W laSdV )> I 5 ' 15 ) 

one sees that it is proportional to this gradient (exterior derivative). 

The Killing- Yano tensor [15] for this spacetime is naturally defined by using the 
principal null directions of the Riemann tensor. The principal Newman-Penrose tetrad 
is 

r 2 + a 2 + I 2 a 
I = -r d t + d r + —d^ , 

' [(r 2 + a 2 + l 2 )8 t + Ad r ~ ad^} , 



2[r 2 + (l + ac) 



' [Zdt - idg + -d$] . (5.16) 



V2[(l + ac) -ir)'s 
The Killing- Yano 2-form is then [15] (their Eq. (35.64)) 

/ = —(/ + ac)[l A n] + ir[m A fh] . (5.17) 
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By taking its self-dual combination K. = / + i*f, one finds 
K, = i[r + i(£ + ac)][dt A dr - iasdt A d6 + xdr A dcp - is(r 2 + a 2 + l 2 )d6 A d<j>] 
= - 4l {lr-^ + *c)}(M + i l)~Q 2 } 

[r + i(£ + ac)}*[r - i(l + ac)} ' 1 ' 

namely it is proportional to the self-dual combination l|3.6(l of the rescaled Papapetrou 
2-form H3.6fl associated with the timelike Killing vector. 

The (symmetric) Cotton- York tensor [28] for the quotient-space geometry can be 
evaluated using either spatial metric since its covariant form is conformally invariant. 
For the rescaled metric (see [10]), this tensor is 

y = -Scurl( (3) Ricci) (5.19) 

and its only nonzero components are 

, _ A cot 9 6sa 2 (e 2 +M 2 -Q 2 )(r-M) 
V -T~M V ' V ~ (A-aVf (5 - 20j 

The condition £ 2 +Ai 2 — Q 2 = makes this tensor zero, corresponding to the conformal 
flatness of either spatial metric. The same is true for the static limit a = 0, where the 
Cotton- York and Simon tensors coincide. 



6. Conclusions 



The gravitoelectric description of the Simon and Simon-Mars tensors which play a 
crucial role in characterizing the Kerr spacetime has been generalized to the class of 
electrovacuum spacetimes where the corresponding Kerr-Newman spacetime exhibits 
a similar behavior once the Simon tensor is suitably generalized to incorporate the 
electromagnetic field. Again the alignment of the principal null directions of the 
Weyl, Papapetrou and now additionally the electromagnetic field is equivalent to the 
vanishing of the Simon-Mars tensor, which in turn forces the generalized Simon tensor 
to vanish as well. The Scurl and cross-product notation and their properties also help 
illuminate some of the manipulations of exact solution theory whose motivation is not 
so obvious. 
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Appendix A. Simon tensor for stationary vacuum spacetimes 

If h a b is the spatial metric, namely the Riemannian metric on the quotient 3-space 
of a stationary spacetime by its Killing trajectories, then a convenient conformal 
rescaling h a b = M 2 h a b by the absolute value of the norm of the timelike Killing 
vector field M 2 = -£ a £ Q > (where it remains timelike) leads to simpler formulas 
for the spacetime Ricci tensor. Like the original definition of Simon, the formula due 
to Perjes [29] quoted as Eq. (18.70) in [15] relies on this conformal metric to express 
the Simon tensor in terms of the Ernst potential. We denote this tensor by SimonES 
(for the book title Exact Solutions) with a tilde (since its indices should be raised and 
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lowered with h ab ) to distinguish it from our Simon tensor defined slightly differently 
using the original spatial metric 

SimonES b a = [Me£}- 2 f} bcd [£ c V d £ a - h ac h f9 £ [d V g] £ /] , (A.l) 

where £ a = V a £ = V a £ = £ a , $le£ = —M 2 , and all the metric and duality operations 
refer to the conformally rescaled 3-metric. (In this appendix, V refers to the covariant 
derivative associated with the metric h a b-) Apart from an extra factor of 2, Krish [30] 
also uses this definition. 

One easily sees that this tensor is symmetric, and since the second term is 
antisymmetric, it must exactly cancel the antisymmetric part of the first term, so 
only symmetric part of the first term remains and that is just the Scurl operator 
acting on the gradient of the Simon potential apart from the scalar multiplier, which 
is automatically tracefree. Thus the Simon tensor is much more efficiently represented 
by a formula which makes these properties obvious 

SimonES ah = [3ie £]- 2 f, cd {a V \ c \£ b) £ d = -[5ie£]- 2 [V£x V(V£)U 

= [5ie£]- 2 [Scurl(V£<8) W£)] ab , (A.2) 

where the last equality follows from the gradient identity IjB.ljl in the conformal 
geometry. 

Finally one must use the conformal transformation properties of the Scurl 
operator to re-express this in terms of the original spatial metric. For a symmetric 
covariant tensor X ab = X/ a £X/ b £ which does not transform X ab = X ab , one finds 

[ScurllU = [ScurltM- 1 ^)]^ , (A.3) 

so 

SimonES a6 = [$e£]- 2 [Scurl (M _1 V£ <g> V£)] a6 = AM' 1 SimonMars o6 , (A.4) 

where the last equality follows from Eq. I|4.11|l . 

Simon introduced a well-known special transformation w = (1 — £)/(l +£) of the 
Ernst potential which has the property that it only changes the Scurl combination 
above by an overall factor 

Scurl (Vw Vw) = ~(1 + w) 4 Scurl (V£ O V£) . (A.5) 

As discussed in appendix B, this is actually true for any function of the Ernst potential. 
The spatial dual of the original definition of Simon 

C bcd = 4(1 - wwy 2 (V [d wV c] V b w - h b[c U d] ) , U d = h fg V [d wV g] V f w , (A.6) 
can therefore be rewritten in the same way as the above definition 

\fla Cd C bcd = 2(1 - ww)- 2 fj a cd (V [d wV c] V b w - h bc h f 9 V [d wV g] V s w) 

= 2(1 - wwy 2 [Scurl (Vw <g> Vw)] ab 

= 2(1 - ww)~ 2 (l + wf [Scurl (V£ <g> V£)] a b 

= 2(1 - ww)- 2 (l + w) i [Scuil(M- 1 V£ <E> V£)] ab . (A.7) 

The Scurl formula for the Simon tensor expressed in terms of the rescaled 
metric h ab has an obvious immediate consequence for the interpretation of its 
vanishing. Because of the Scurl identity l|B.2|) in the conformally rescaled geometry, 
the differential condition Scurl(V£ ® V£) = is equivalent to the algebraic condition 

V£xV(V£) = 0, (A.8) 
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which in turn expresses the "alignment" of the two tensors V(V£ ) and V£ ® V£, i.e. 

[V(V£)] (TF) oc [V£0V£] (TF) . (A.9) 
This geometrical meaning has been overlooked perhaps because of the unnecessary 
complication in previous formulas. 

For the electrovac case, the conformal transformation formula for a vector 
divergence 

divX = M~ 3 div(MX) (A.10) 
helps show the equivalence of Eq. (|3.12|) with the original Ernst equation (18.38) 
of [15]. Using this formula and substituting for z in terms of £ and $ by solving 
Eq. (|?TTT|) in E q- EUlll ; one obtains 

divS - M~ 2 h ab {£ a + 2$V a $)£& = . (A.ll) 



Appendix B. Scurl and cross product identities 

In this appendix we use the simpler notation V for the spatial covariant derivative and 
assume all fields are spatial and stationary. A straightforward computation verifies 
the identity 

Scurl(X X) = ^ [SYM(A ® curlX)] (TF) - [X x SYM(VX)] (B.l) 

for any vector field X. When the vector field is a gradient X — V/, then 
curlX = curl(V/) = and the covariant derivative VI = V(V/) = SYM(V(V/)) is 
symmetric, so this reduces to the useful identity (see Eq (3.19) of [10]) 

Scurl(X <g> X) = -X x (VI) if X = V/ . (B.2) 
Notice that for any function F(f) of /, then WF(f) = F'(f)V f and one has the 
property that the rate of change factor passes through the Scurl operation without 
extra terms, expanding the Scurl by the product rule for a scalar and a tensor 
Scurl(VF ® VF) = Scurl(F' 2 V/ <g> V/) 

= (F' 2 )'Wf x (V/ «) V/) + F' 2 Scurl(V/ ® V/) = F' 2 Scurl(V/ <E> V/) (B.3) 
because of the identity Y x (Y (g) Y) = valid for any vector field Y. 

When performing actual calculations in this approach, one needs to refer to a 
number of Scurl and cross product identities involving vector fields and symmetric 
tensors where the symmetric tensors themselve arise from vector fields. For arbitrary 
spatial vector fields X, Y and Z one finds by straightforward evaluation 

Scurl[SYM(A ® Y)) = - -[X x SYM(VF) + Y x SYM(VA)] + 

+ - SYM[(X <g> curlY + Y <g> curlX)]( TF ) , (B.4) 

Scurl[SYM(VA)] = - [X x ^Ricci] + - SYM[V(curl A)]( TF ) , (B.5) 
SYM[V(X x Y)] {TF) = X x SYM(W) -Fx SYM(VA) + 

+ i SYM[A curlF -78 curbY] (TF) , (B.6) 

SYM[X (Y x X)} =Y x[X®X] , (B.7) 

ixSYM[i®y] =-hrx[x®x], (b.8) 

XxSYM[YxZ] = - SYM[y ® {X x Z) + Z ® {X x Y)] . (B.9) 
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